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A numerical study of Orientifold Planar Equivalence is performed in SU(N) Yang-Mills theories 
for N = 2,3,4,6. Quenched meson masses are extracted in the antisymmetric, symmetric and 
adjoint representations for the pseudoscalar and vector channels. An extrapolation of the vector 
mass as a function of the pseudoscalar mass to the large-rV limit shows that the numerical results 
agree within errors for the three theories, as predicted by Orientifold Planar Equivalence. As a 
byproduct of the extrapolation, the size of the corrections up to 0(1 /N^) are evaluated. A crucial 
prerequisite for the extrapolation is the determination of an analytical relationship between the 
corrections in the symmetric and in the antisymmetric representations, order by order in a 1 /N 
expansion. 
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1. Introduction 

Orientifold Planar Equivalence [jl|] is a powerful analytical tool that establishes the equality of 
certain observables (among which, meson masses) in two classes of gauge theories at N = °°, N 
being the number of colours. This equivalence is of particular interest when the two gauge theory 
families it relates are S\J(N) with Nf (anti)symmetric Dirac fermion flavours and SU(A0 with Nf 
Majorana flavours in the adjoint representation. In fact, this special case enables us to transcribe 
SUSY results to QCD [Q], provided that the latter is close (in the sense of the 1 /N expansion) to its 
large-Af limit. Some issues related to the equivalence, like e.g. the size of the finite N corrections, 
mandate an ab-initio calculation. A full dynamical calculation is expensive from a computational 
point of view. In fact, the cost of the computation is brought almost entirely by the inversion of the 
fermionic matrix, which for theories with two-index representation fermions requires 0(N 4 ) oper- 
ations. Before undertaking a full dynamical study, it is convenient to consider the quenched case 
as a prototype example. Although for theories in two-index representations the quenched theory 
and the dynamical theory do not coincide in the large-Af limit, Orientifold Planar Equivalence still 
holds in the quenched case, which then becomes a useful toy model to study techniques that could 
be used in the dynamical investigation. In this spirit, a first paper appeared in which the chiral 
condensate was shown to coincide in the large-N limit of the adjoint, symmetric and antisymmet- 
ric representations [Q]. A crucial point in that work was the separation of the corrections in even 
and odd powers of 1 /N by combining the numerical results in the symmetric and antisymmetric 
representations. As a result, a precise extrapolation to the large-Af limit could be performed and 
Orientifold Planar Equivalence was verified to a high degree of accuracy. Moreover, this technique 
allowed a computation of corrections up to 0(1 /N 3 ) in the (anti)symmetric representation while 
still using a reasonably small value of N. In this work, we will perform a similar study for the 
quenched meson spectrum. First, we will show that there is a simple order by order relationship 
between the coefficients in a large-Af expansion of amplitudes in correlation functions and of spec- 
tral masses in the symmetric and antisymmetric representations. Then, we shall use a simple chiral 
ansatz to establish a relationship between the mass of the vector meson my and the mass of the 
pseudoscalar meson m PS at finite N, for N ranging from two to six, in the adjoint, in the symmetric 
and antisymmetric representations. Finally, using the analytical relationships between the correc- 
tions, the chiral ansatz will be extrapolated to the large-Af limit and Orientifold Planar Equivalence 
proven to hold. The results we will present in this work have already appeared in [Qj, to which we 
refer for a more detailed discussion. 

2. Meson masses for two-index representation fermions in the 1 /N expansion 

We extract meson masses from the large distance exponential decay of the meson correlation 
functions. For the fermion sources we use the Wilson formulation, in terms of which the Dirac 
operator is given by 

D^= (m + tySxy-Kxy, (2.1) 

with m the bare quark mass and 

Kxy = ~\ [{r-^)R[U^x)} 5,^ + {r + J^R [l/fa)] S, X _ A ] . (2.2) 
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In the previous formulae, x and y are site indices and the the Dirac T matrices in Euclidean space. 
R [Ufi(xj\ is the link variable stemming from x in the direction jx expressed in the representation R. 
In our case, R can be the adjoint (Adj), the two-index symmetric (S) or the two-index antisymmetric 
(AS) representation. For the representation R, meson correlators take the form 

cfV 2 fry) = r * ( tr « (^r^v 1 ^) ) yM , (2.3) 

where tr# indicates the trace over the color indices in the representation R and T\ and T2 are 
combinations of T matrices associated with the quantum numbers J PC of the mesons. With the 
choice r« = 1 for R = S/AS and r« = 1/2 for R = Adj, in the large-A^ limit we have the equality of 
correlators in these three representations. The factor r# = 1/2 in the adjoint representation can be 
understood considering that in our formulation we are always taking Dirac flavours in the adjoint, 
and one Dirac flavour is equivalent to two Majorana flavours. 

The formal proof of the equality of the correlators can be sketched as follows. We start ex- 
panding the correlators in Wilson loops: 

-LqK ra (x,y) = £ £ ay<tr*WV), (2.4) 

'tf being a closed curve that contains both x and y. The coefficient a%> in (14) does not depend on 
the representation. Writing explicitly the traces in two-index representations in terms of the trace 
of the fundamental representation yields 



2\ 



^ = 2^* 

WD(x,y) 



(2.6) 



Note that the factors of 1/2 in the two previous formulae have a different origin: for the S/AS, it 
comes when expressing the trace in the two-index representations in terms of the group element 
in the fundamental representation, while for the adjoint is the factor tr. At large N, we can use 
factorisation and neglect the subleading terms, obtaining 

'4'?fry) = ~ I or^P, (2-7) 

1 Adj ( , 1 v n (trWV)(trW$) 



Finally, the equality 



follows from charge conjugation invariance, which dictates (trW^.) = (trW<g). This last step in 
Yang-Mills theories does not present the conceptual issues discussed in \K H] for the dynamical 
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case. The equality of the meson spectra in the three theories follows from the equality of the 
correlators. 

Correlators of mesons in the adjoint representation are expected to be expandable in a power 
series of \/N 2 . As for the other two representations, since in the symmetric and antisymmetric 
representations the static sources introduce 1 /N effects, there is no reason not to expect 1 /N cor- 
rections in the meson mass spectrum. Assuming that there is no accidental degeneracy at large-Af 
in the theory with fermions in the antisymmetric and in the symmetric 1 , the masses Mr and the 
amplitudes Ar in correlators of fermion bilinears in these two representations can be shown to be 
related by 

A } s (N)=A j AS (-N) and m J s (N) = m J AS (-N) (2.10) 

for states j related by Orientifold Planar Equivalence. Order by order in a l/N expansion, this 
implies that even-power terms in 1 /N in the symmetric and antisymmeuic representation have the 
same coefficient, while odd-power terms have opposite coefficients. This means that the combina- 
tions 

M j = (m{ + m j AS ) /2 and \x j =n(m } s - m j AS ) /2 (2.1 1) 

have an expansion in powers of l/N 2 . 



3. Numerical results 



We generated gauge configurations for values of N ranging from two to six using the standard 
Wilson action. Both the generation of the configurations and the analysis of the correlators in 
two-index representations were performed using the purposely developed HiRep code [0|. The 
technical details of the simulations have been given in [H|. Here we remark that we simulate at 
the same j8 values as in Ref. [Q], where the couplings were chosen so that the gauge theories at 
various N are at a common value of the lattice spacing a = l/(5T c ), T c being the deconfinement 
temperature of the pure gauge system [|8j|. 

After extracting my and mps, we have determined the chiral limit value of my (denoted as my) 
using the chiral ansatz 

my(mps)=my+Bm PS . (3.1) 

Our results and the corresponding fits are displayed in Figs. [T], |2| and ^respectively for the adjoint, 
symmetric and antisymmetric representations. We note that there is a well-defined hierarchy for my 
at fixed mp$, with the symmetric mass being the highest and the antisymmetric being the lowest. 
At fixed representation, the vector mass in the antisymmetric has the largest variation with N, while 
there is little or no variation with N for the adjoint vector masses. 

The large-Af extrapolation has been performed for the quantity rn yAd . and for the quantities 
My and fly (see Eqs. ( 2.1 ![ )) using in all cases only the leading term and the 1 /N 2 correction. The 



'in Quantum Field Theory, degeneracies are due to a symmetry, and there is no symmetry that could generate 
degeneracies in these theories. 
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Figure 1: my as a function of m 2 ps for the adjoint representation at the values of N shown. The solid lines 
are chiral fits to the data. 
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Figure 2: my as a function of rrip S for the symmetric representation at the values of N shown. The solid 
lines are chiral fits to the data. 



quality of the fits is displayed in Fig. |]. The values of m,y S and «Jy AS as a function of \/N are 
obtained inverting Eqs. ( [2.1 1[ ). Our results are: 

Adj= 0-736(17) (3.2) 
m ^„, 2 3( 27) + ^_i^ + ^, 

4iS .o, 2 3 P7) -!^-i^-^). 

The equality of the large-/V results for the three representations is a prediction of Orientifold Planar 
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Figure 3: my as a function of m^ s for the antisymmetric representation at the values of N shown. The solid 
lines are chiral fits to the data. 
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Figure 4: Large-Af extrapolation of the adjoint mass and of the combinations M and fi, all in the chiral limit. 



Equivalence. We remark that through the separations of corrections in even and odd powers of 1 /N 
in the symmetric and antisymmetric representation, we were able to get an handle on the size of 
the corrections up to 0(1 /N 3 ). Numerically, Orientifold Planar Equivalence can be verified for any 
value of mps in the chiral region 

4. Conclusions 

Our work should be regarded as an illustration of a possible strategy for a numerical proof 
of Orientifold Planar Equivalence. We have shown that it is possible to combine the data in the 
symmetric and antisymmetric representation for extracting precise values in the large-Af limit from 
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data at moderately small values of N. The next step in this investigation is the simulation of the 
dynamical theories. 
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